Fermi-liquid effects in the gapless state of marginally thin superconducting films 
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We present low temperature tunneling density-of-states measurements in Al films in high parallel 
magnetic fields. The thickness range of the films, t — 6 — 9 nm, was chosen so that the orbital and 
Zeeman contributions to their parallel critical fields were comparable. In this quasi-spin paramag- 
netically limited configuration, the field produces a significant suppression of the gap, and at high 
fields the gapless state is reached. By comparing measured and calculated tunneling spectra we are 
able to extract the value of the antisymmetric Fermi-liquid parameter G'^ and thereby deduce the 
quasiparticle density dependence of the effective parameter Gj!ff across the gapless state. 

PACS numbers: 74.78.Db, 73.40.Jn, 74.25.-q 



I. INTRODUCTION 

One of the most successful constructs in condensed 
matter physics is Landau's Fermi-hquid (FL) description 
of the low temperature properties of interacting electron 
systemsiii^ The theory provides for the salient thermody- 
namic and transport properties of a profoundly complex 
many-body system in terms of a low-density gas of quasi- 
particle excitations. The averaged effects of electron- 
electron and electron-phonon interactions are incorpo- 
rated into a number of FL parameters, which are as- 
sociated with the renormalization of fundamental sys- 
tem properties such as the effective mass and the spin 
susceptibility— In particular, the so-called antisymmetric 
FL parameter G° affects the spin response of the system 
and is related to the ratio of the spin susceptibility den- 
sity of states N{x) to the heat-capacity density of states 
7V(7) by G° = N{j)/N{x) - l^^ In this paper we report 
low-temperature measurements of the parameter G'^ in 
superconducting disordered Al films. 

Although is a fundamentally important parameter 
in the many-body description of metals, there have been 
very few measurements of it reported in the literature. 
This is due, in part, to the fact that it is difficult to 
measure it directly in bulk systems. Consequently, little 
is known about its dependence on disorder, quasiparti- 
cle density, and/or dimensionality. Measurements of 
have been extracted from high-field tunneling density of 
states (TDOS) and critical-field studies of low atomic 
mass superconducting films. For example, estimates of 
normal-state value of G° have been reported via tunnel- 
ing studies of the Zeeman splitting of the BCS density of 
states in superconducting Al films, G° ^ 0.3 — 0.4,'*'^ in 
the intermediate to high-temperature regime where the 
superconducting order parameter is partly suppressed 
by thermal fluctuations. Low-temperature TDOS mea- 
surements in amorphous Ga films, which are strong- 
coupling superconductors, give a somewhat larger value 
of G° ~ 0.81j^ Alternatively, G° can be extracted from 
parallel critical field measurements. This method gives 



G" ^ 0.23 in TiN films^^i The only direct measurements 
of G" have been obtained in the normal state, via the field 
dependence of the pairing resonance* in Al, G° ^ 0.17;^ 
and Be films, G° ^ 0.21.^*^ In all these systems the spin- 
orbit scattering is quite low; thus, spin remains a good 
quantum number. 

In the experiments described above the films were suf- 
ficiently thin so as to suppress the Meissner currents; 
therefore the field response was purely due to the elec- 
tron Zeeman splitting (see e.g. Ref. \lu and Sec.|TT|. By 
contrast, in this work we consider thicker films; this en- 
ables us to explore the gapless state and determine G° 
from low temperature TDOS measurements. In the nor- 
mal state, the internal field Hi felt by the electron spins is 
given by Hi = H + ifex, where H is the applied field and 
the exchange field is ffox = -HG°/{1 + G°). Deep in the 
superconducting phase, i.e, at low temperature and low 
field, the number of quasiparticles is small due to the fact 
that they have been consumed by the formation of the 
superconducting condensate. In this limit, the exchange 
effects are greatly suppressed and the effective Landau 
parameter [see Eq. ^] is small compared to its normal 
state value, G^g ^ 0. However, as one approaches the 
critical field Hc\\, the quasiparticle density increases and 
G°g — > G°. Below we show that the main features of the 
TDOS in the gapless superconducting state are strongly 
influenced by the rise of these exchange interactions. In 
addition to its fundamental interest, a thorough under- 
standing of exchange effects in superconducting Al films 
is important for possible applications in mesoscopic hy- 
brid structures aimed at the control of current spin polar- 
ization via the Zeeman splitting in the superconducting 
elements. 

This paper is organized as follows: in Sec. |TT] we give 
a brief overview of the phase transition in thin films and 
introduce the parameters necessary for the description 
of its properties. In Sec. IIIII we give some details of the 
samples' preparation, and in Sec. llVl we discuss our main 
results. 
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II. PHASE TRANSITION IN MARGINALLY 
THIN FILMS 

In bulk superconductors the critical-field transition is 
completely dominated by the orbital response of the con- 
duction electronsp^ It is possible, however, to inhibit the 
orbital response by applying a magnetic field in the plane 
of a superconducting film whose thickness t is much less 
than the superconducting coherence length ^ and whose 
electron diffusivity is lowJ^^^^^ Under these conditions 
the phase transition to the normal state is mediated by 
the spin polarization of the electrons, and at the criti- 
cal field the electron Zeeman splitting is of the order of 
the superconducting gap energy. At zero temperature, 
in particular, a first-order transition from the supercon- 
ducting state to the paramagnetic normal state occurs at 
the Clogston-Chandrasekhar critical fiel d^^'^^ 



ttCC _ 



(1) 



where Aq is the zero-field, zero-temperature gap energy, 
and is the Bohr magneton. This so-called "spin- 
paramagnetic" (S-P) transition is realized in low atomic 
mass superconductors such as Al, Be, and TiN, and it is 
accompanied by very peculiar dynamics i^ i-'^^'^^i^^'^" Note 
that for positive G° the parallel critical field exceeds the 
non-interacting Clogston-Chandrasekhar limit. 

In the thin- film limit and at low temperatures, the par- 
allel magnetic field has little effect on the order parame- 
ter up to the transition, and the Clogston-Chandrasekhar 
critical field [Eq. ([T])] can be more than an order of mag- 
nitude higher than its orbitally-mediated counterpart. 
However, if the thin film condition is marginally relaxed 
then orbital effects will generally lower the critical field. 
A measure of the relative weight of the orbital response 
compared to the spin polarization in determining the par- 
allel critical field is given by the dimcnsionless orbital 
pair-breaking parameterli 



(2) 



where D is the electron diffusivity, t is the film thickness, 
e is the electron charge, and £ is the mean-free path for 
impurity scattering. The function f{i/t) describes the 
crossover between local electrodynamic response of the 
film for £ ^ t, where f{(/t) 1, and the non-local one 
in the opposite limit £ > t, where f{e/t) ~ ^t/ii^ For 
marginally thin films used in this study, c ^ 1. 

Beside the film's thickness (through the parameter c) 
and exchange effects, spin-orbit scattering also affects 
the value of the parallel critical field. Therefore, fol- 
lowing Ref. [Ill, we introduce a third dimensionless pa- 
rameter b = h/SAoTso, where 1/tso is the spin-orbit 
scattering rate. As remarked in Sec. [H in the materi- 
als under consideration spin-orbit scattering is a small 
effect, 5 <C 1. The three parameters c, G°, and 6, 



can have competing effects on the value of the paral- 
lel critical field. Moreover, they also affect the position 
of the tricritical point separating the low-temperature 
first-order phase transition from the higher-temperature 
second-order one. In the limit G'^jC ^ the tricritical 
temperature is Ttri — 0.56Tc, where Tc is the critical 
temperature.^^ More generally, the tricritical tempera- 
ture and other properties of a superconducting film, such 
as the TDOS, can be calculated by solving the Usadel 
equations for the semiclassical Green's functions together 
with the self-consistent equations for the order parame- 
ter S and the internal magnetic field Hi] these equations 
can be found in Ref. (see also Ref. iTj for an alterna- 
tive parametrization). All three dimensionless parame- 
ters enter into these mean field equations, while Aq (or 
Tc ~ O.STAo/fcs) fixes the overall energy scale. 

In this work we present low temperature measurements 
of TDOS in Al films whose thicknesses are about two to 
three times greater than the t ~ 3 nm typically used in 
spin-paramagnetic studiesi^ Though the thickness of the 
films remains much smaller than the superconducting co- 
herence length, t ^ ^, the orbital response to the applied 
field is no longer negligible. Indeed, the thickness range 
is chosen to assure that the orbital and Zeeman contri- 
butions to the critical field are comparable (this could 
also be accomplished by rotating the films out of parallel 
orientation, but a finite perpendicular field component 
would introduce vorticity^). Interestingly, the critical 
field transition can become reentrant^^ under these con- 
ditions. The reentrance is associated with a high-entropy 
gapless state near the critical field, in which the supercon- 
ducting state has higher entropy that the normal statej^ 
The primary purpose of the present work is to determine 
through TDOS measurements across the gapless re- 
gion, where there is a monotonic increase in quasiparticle 
density as one approaches the second-order phase transi- 
tion to the normal state. 



III. SAMPLE PREPARATION 

Aluminum films were fabricated by e-beam deposition 
of 99.999% Al onto fire polished glass microscope slides 
held at 84 K. The depositions were made in a typical vac- 
uum P < 3 X 10^^ Torr at a rate of ^ 0.2 nm/s. Films 
with thicknesses ranging from 6 to 9 nm had normal- 
state sheet resistances that ranged from i? = 10 to 20 il 
at 100 mK. After deposition, the films were exposed to 
the atmosphere for 0.5-4 h in order to allow a thin native 
oxide layer to form. On top of the oxide, serving as the 
tunneling barrier, a 14-nm-thick Al counterelectrode was 
deposited. Due to its relatively large thickness, the coun- 
terelectrode had a parallel critical field of ~ 1.1 T, which 
is to be compared with Hc\\ ~ 3 T for the films. The 
junction resistance ranged from 10 to 20 kfi, depend- 
ing on exposure time and other factors, for a junction 
area of about 1x1 mm^. Only junctions with resistances 
much greater than that of the films were used. Magnetic 
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FIG. 1: Upper panel: Calculated field dependencies of the su- 
perconducting order parameter S, gap A, and ratio between 
internal and applied fields Hi/H for a 2.5-nm-thick Al film 
(Tc = 2.7 K) in parallel field at 77 mK. The arrow indi- 
cates the theoretical supercooling field. Note that the gap 
is finite up to the superheating field. The upper inset is the 
corresponding phase diagram where S is the superconducting 
phase, N the normal state, and SM is the state memory re- 
gion. Lower panel: Field dependencies of the same quantities 
for a 7-nm-thick Al film {Tc = 2.1 K) at 60 mK. The verti- 
cal dot-dashed line separates the gapped and gapless regions. 
Note in the latter the fast drop of the order parameter with 
increasing field. The lower inset is the corresponding phase 
diagram, with a maximum in the second-order critical-field 
curve. The two points in the phase diagram represent the 
temperature and fields at which tunneling spectra were taken 
in the gapless region. 



fields of up to 9 T were applied using a superconducting 
solenoid. An in situ mechanical rotator was employed to 
align the film surface paraUel to the field with a preci- 
sion of ~ 0.1°. Measurements of resistance and tunnehng 
were carried out in an Oxford dilution refrigerator using 
a standard ac four-probe technique. 



IV. RESULTS AND DISCUSSION 

Shown in the upper panel in Fig. [T] are the calculated 
low-temperature field dependencies of the order param- 



eter 5, gap energy A, and internal field Hi for a 2.5- 
nm-thick Al film having = 2.7 K, AJe = 0.41 mV, 
^ ^ 15 nm, and sheet resistance R ^ 1 kSl. At this thick- 
ness and resistivity orbital effects are negligible, since we 
estimate c ~ 0.02 ^ 1. Accordingly, 5 and Hi/H are rel- 
atively insensitive to the applied field up to the first-order 
parallel critical field located approximately half way be- 
tween the supercooling field (arrow in figure) and the 
superheating one (used for normalization in the upper 
panel). In contrast the gap A, while remaining finite, 
decreases linearly as expected because even a small spin- 
orbit scattering mixes the densities of states of opposite 
spins, which in turn are shifted in opposite directions 
by the Zeeman field. The curves are computed by nu- 
merically solving the mean-field equations mentioned in 
Sec. ini for concreteness, we use for the spin-orbit scat- 
tering parameter b — 0.052, and for the Fermi-liquid pa- 
rameter = 0.155, as these values are the one wc ex- 
tract from the measurements in thicker films. Shown in 
the inset of the upper panel of Fig. [1] is the measured 
parallel critical field of this film as a function of temper- 
ature. This plot represents the classic S-P phase diagram 
in which a high-temperature line of second-order phase 
transitions crosses over into a line of first-order transi- 
tions at the tricritical point Ttri ^ O.STcJ^ 

In the lower panel of Fig. [1] we show the corresponding 
behavior of (5, A, and Hi at 60 mK for a 7-nm-thick Al 
film having ~ 2.1 K, A^/e ~ 0.32 mV, ^ - 60 nm^ 
and R = 16 fi. Not only does this thicker film have a 
lower critical field due to its finite orbital response, but 
also the structure of its phase diagram (lower inset) is 
dramatically altered from that of the classic S-P dia- 
gram. Indeed, due to the local maximum in the parallel 
critical field near T/Tc — 0.45, the critical-field behav- 
ior is reversibly reentrant. The most obvious difference 
between the field dependencies shown in the upper and 
lower panels of Fig.[l]is that the gap energy in the thicker 
film goes to zero before the critical field is reached. In 
fact, the region to the right of the vertical dot-dashed 
line represents a gapless superconducting state; clearly, 
in this region rapid changes in 5 and Hi take place. 

Deep in the superconducting phase, where the super- 
conducting order parameter is well established, exchange 
effects are negligible due to the fact that the quasiparticle 
density is low and the e-e interaction effects parametrized 
by are preempted by the formation of the conden- 
sate. Interestingly, however, it is clear in both panels 
of Fig. [T]that Hi/H never reaches unity in the limit of 
small applied field. We believe that this is a spin-orbit 
effect, and that, in the zero-temperature zero- field limit 
Hi/ H = \/ \l + G^a{h)\ , where a{h) is a positive func- 
tion of the spin-orbit parameter h. a{b) goes linearly to 
zero as 6 — > and saturates to 1 in the limit b —^ oo. 
Once the field is increased beyond the onset of the gap- 
less state, the quasiparticle density rapidly grows until 
it reaches its normal-state value at iJ^i as we discuss 
in more detail below. Similarly, exchange effects, as re- 
flected in the internal field, also begin to "turn on" upon 
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FIG. 2: Tunneling conductance normalized by the normal- 
state conductance G„ as a function of normalized bias voltage. 
The solid lines are experimental spectra of a 7-nm-thick Al 
film taken at 60 mK in a parallel field of 2.20 (upper panel, 
gapped) and a 3.02 T (lower panel, gapless). The dashed 
curves are the calculated spectra - see the text for more de- 
tails. The arrows point to the primary peak and dip features 
of the spectra. The insets show spectra for an 8-nm-thick film 
taken at 1.40 (upper inset) and 2.42 T (lower one); the axes 
cover the same ranges as the respective main panels. The 
dashed lines axe the theoretical curves hardly distinguishable 
from the continuous experimental lines. 



entering the gapless state. In contrast, in the purely S- 
P limit, there is no gapless state, and the internal field 
jumps discontinuously to its normal-state value at the 
first-order critical-field transition (see the upper panel of 
Fig.[T]). 

In Fig. [2] we compare the normalized tunneling con- 
ductances obtained in the gapped and gapless states of 
the 7 nm Al film at 60 mK. The upper panel was taken 
at _ff = 2.2 T corresponding to the gapped phase. The 
curve clearly displays a Zeeman splitting of the BCS 
density-of-states peaks.— The data in the lower panel 
of Fig. [H were measured at 3.02 T. Note that the gap is 
completely suppressed and that a large number of states 
exist at the Fermi energy {i.e., V=0). The dashed lines 
in Fig. [2] are the theoretical curves calculated with pa- 
rameters c = 0.79, G° = 0.155, and b = 0.052.25 We also 



find good agreement between theory and experiment for 
TDOS measurements (see insets of Fig. [2]) in an 8-nm- 
thick film with c ~ 1.23; this value is in reasonable agree- 
ment with the scaling c (x oi Eq. ^ in the non-local 
limit relevant to these samples. jj^ the lower panel, an 
additional broadening,— F = O.OI6A0 ~ 5 x 10^® eV, is 
used in calculating the theoretical curve; i.e., the rela- 
tion between density of states v and Green's function Q 
is taken to be J'(e) oc Im 0{e + iT). A broadening of sim- 
ilar magnitude (F = O.OIAq) is used for the calculated 
curves in both insets. This finite broadening is larger 
than values reported in the literature)^ but at these fields 
it could be caused by a small misalignment of the sam- 
ple. Indeed, the Cooper pair cyclotron frequency asso- 
ciated with a field misalignment of angle 9 with respect 
to the sample plane is hnc± — 4eDH sm9; a,t 6 = 0.1° 
and H = 3 T, and using the estimate D ~ 25 cm^/s, 
we obtain hQc± — 5 x 10~^eV > F. We also note that 
small variations in the parameters c and F can compen- 
sate each other. On one hand this seems to support the 
idea that in the present case F could be caused by an 
orbital effect of the field as c is; on the other hand this 
interdependence makes it difficult to extract F from the 
daXaM- and a more thorough study of the broadening, 
including low-field measurements and/or intentional tilt- 
ing, is necessary to verify this interpretation. 

Although the two tunneling spectra in Fig. [5] are very 
different in structure, the evolution of the spectra as one 
moves from the gapped to the gapless state is contin- 
uous. The salient peak and dip features of each spec- 
tra are denoted by the arrows in Fig. [2] (we point out 
that the small broadening F, introduced above, affects 
the amplitude but not the position of these features in 
the theoretical curves). The behavior of these features as 
one crosses over into the gapless state is shown in Fig. |31 
The main panel depicts the low-temperature field depen- 
dence of the outermost coherence peaks along with data 
points corresponding to tunneling measurements. Inter- 
estingly, the position of the peaks is not monotonic in 
field, and the peaks move to lower energy in the gapless 
state. The dotted line is simply the Zeeman splitting, 
which represents the expected peak position if there were 
no orbital depairing and no exchange effects. We have 
verified numerically that, in contrast to the results for 
the marginally thin films, the peaks follow the bare Zee- 
man splitting with good approximation for the thinner 
film with c = 0.02, in agreement with low-temperature 
measurements in thinner films The dashed curve is 
the expected peak position with orbital depairing but 
no exchange. This curve misses the data points badly 
while agreeing with the solid curve at low field indicating 
that G*^ plays a significant role in determining the tun- 
neling density of states in the gapless state. Moreover, 
with — the low-temperature phase transition is pre- 
dicted to be first-order, in contrast to the experimental 
finding of a second-order transition. The inset of Fig. [3] 
shows the corresponding field dependence of the spectra 
dips. As is the case with the peaks, the evolution of the 
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FIG. 3: Position of the peaks in Fig.[2]as a function of parallel 
field. The solid line represents the theoretical peak positions 
with Fermi-liquid effects included. Note that in the gapless 
state (right of the dot-dashed vertical line), the peaks move to 
lower energy as the field increases; this behavior is in agree- 
ment with the experimental measurements (circles with error 
bars) . The dashed curve shows the theoretical peak positions 
in the case of no exchange, i.e., = 0. The dotted line 
represents the Zeeman energy, which is the expected peak 
position with no exchange and no orbital pair breaking. In- 
set: Dip positions as a function of field. The solid line is the 
theoretical position including exchange, and the dashed line 
with exchange excluded. Triangles are data points. 



dips' position contrasts with that of the thin-film case. 
In the latter, the dips follow the peaks to higher energies. 
In Fig. [3l the dips move to lower energy and the peaks 
broaden with increasing field. This broadening could po- 
tentially limit the degree of spin polarization that could 
be achieved in hybrid structures 

Building on the above analysis we can now extract an 
"effective" PL exchange parameter G'^g as a function of 
the quasiparticle density n„. In Fig. |3] we plot such a 
function, where we assume a two-fiuid model in order 
to estimate n„. In this phenomenological approach, the 
superconducting carrier density rig is related to the pen- 
etration depth A by 



cx A 



(3) 



and the normal electron density is n„ = 1 — Ug. In disor- 
dered thin superconducting films the parallel field pene- 
tration depth is proportional (up to a weakly tempera- 
ture dependent coefficient) to the penetration depth Al 
of a clean superconductor in the local limit [see, e.g., Eq. 
(3.136) of Ref. [ij. The latter depends on temperature 
approximately as 



Al(7) oc 
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FIG. 4: Effective FL parameter [Eq. ((6)1] as a function 
of quasiparticle density n„ in the superconducting phase of a 
marginally thin Al film. The inset is the calculated quasipar- 
ticle density as a function of the normalized parallel magnetic 
field for a 7-nm-thick Al film at 60 mK. The vertical dot- 
dashed line divides the gapped and gapless states. Note that 
in the latter region, the quasiparticle density quickly rises 
with the applied field. 



where 7 — T/Tc. Near the critical temperature, the tem- 
perature dependence of the order parameter is (^(7) / Aq c>c 
— 7. Inverting this relationship to express 7 in terms 
of (5, substituting the result into Eq. (|4]), and using the 
definition of n„, we arrive at: 



1 - 



^0 



(5) 



From this relation we evaluate the field dependence of the 
normal electron density through the calculated field de- 
pendence of the order parameter (at low temperature the 
proportionality constant can be taken to be unity up to 
exponentially small corrections); the result is presented 
in the inset of Fig. H) We define the effective Fermi- liquid 
constant by 



Gcff - H/Hi-\ 



(6) 



We can obtain this quantity from the calculated ratio 
Hi/H shown in Fig. [T] Note that G^g is approximately 
linear in the quasiparticle density, thus showing the grad- 
ual rise of exchange effects going from the superconduct- 
ing to the normal phase. This microscopic calculation 
fully supports the qualitative explanation given in Ref.i 
for different g-factors in the superconducting and normal 
states, which correspond to G\ = and G\ ~ 0.17, 
respectively. 

In summary, we have presented measurements of the 
tunneling density of states of marginally thin Al films 
in parallel magnetic field. The marginality condition 
c 1, with the dimensionless parameter c defined in 
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Eq. ([2]), gives access to a gapless superconducting state. 
The inclusion of Fermi-liquid effects, parametrized by the 
quantity G*^ ^ 0.16, is necessary to explain the main 
features (peaks and dips) of the tunneling spectra, see 
Figs. [2] and [3l Although these interaction effects are sup- 
pressed deep in the superconducting phase, they become 
more and more important as one moves through the gap- 
less state to the parallel critical-field transition. In the 
gapless region the effective Fermi-liquid parameter G'^g 
grows quasi-linearly with the quasiparticle density. 
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with (po the flux quantum). 
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For example, no broadening is used in the upper panel 
curve - including it would slightly suppress the peaks' mag- 
nitude. Due to their interdependence, however, different 
values for c and F 7^ could also satisfactorily fit the data. 



